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Abstract. A change of variables is introduced to reduce certain nonlinear stochas- 
tic evolution equations with multiplicative noise to the corresponding deterministic 
equation. The result is then used to investigate a stochastic porous medium equa- 
tion. 
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1. Introduction 

Let U = U{t, x) be a solution of the porous medium equation 

(1.1) Ut = AU"", t>0, X eR, m> 1, 

and let X = X{t) be a semi-martingale. Consider a stochastic version of (11.11) . 

(1.2) du = Au"'dt + udX{t); 

the use of multiplicative noise preserves the positivity of the solution. The first main 
result of the paper is as follows: 

Theorem 1.1. There is a one-to-one correspondence between the solutions of the 
deterministic and stochastic porous medium equations, given by 

(1.3) u{t,x) = h{t)U{H{t),x), 



where 



h{t) = l+ [ h{s)dX{s), H{t)= ! h"'-\s)ds. 
Jo Jo 



Thus, many of the results known for the deterministic equation (11.11) . as summarized, 
for example, in Aronson [I], have a clear counterpart in the stochastic case, and the 
objective of the paper is to the study this correspondence. To keep the presentation 
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from becoming unnecessarily abstract, a particular semi- martingale X is considered: 
X(t) = /g + f{s)dw{s), where w is a standard Brownian motion. The 

results derived in this manner can serve as a benchmark for further investigation of 
the stochastic porous medium equation, in particular, driven by space-time noise. 

Here is a typical consequence of combining (11.31) with the known facts about the 
deterministic porous medium equation. 

Theorem 1.2. Letuo{x) be a non-negative continuous function with compact support 
in M^. Then there exists a unique continuous non-negative random field u = u{t,x) 
such that, for every smooth compactly supported function (f, and every t > 0, the 
equality 

u{t, x)(p{x) dx = / uo{x)(p{x) dx + / / u"^{t^ x)A(p{x) dx 

jRd Jo JR'i 

-\- j / u{t, x)(p{x) dw{t) 
Jo Jr'' 

holds with probability one (in other words, u is the solution of du = Au'^dt-\-udw{t) .) 
In addition, 

(1) The function u is Holder continuous in t, x on every set of the form [T, +oo) x 
R'^, T > 0; 

(2) The mean mass is preserved: 

E / u{t,x)dx= / u{0,x)dx; 
Jr^ Jr^ 

(3) The support of u is uniformly bounded: there exists a random variable rj such 
that, with probability one, < rj < oo and u(t,x) = for \x\ > r] and all 
t > 0. 

(4) For every x G M, \imt-,oo u{t , x) = with probability one. 

Section [2] presents the general result about the change of variables (11.31) for a large 
class of nonlinear equations. Section [3] discusses the basic questions about existence, 
uniqueness, and regularity of the solution for the stochastic porous medium equation, 
and Section H] investigates the long-time behavior of the solution. Because of the 
random time change, the long-time behaviors in the deterministic and stochastic 
cases can be different. 

The following are some comments about the physical origin and relevance of equa- 
tion (11.21) . as well as connections with recent work on the subject. For a variety of 
distributed system, there are two functions that describe the state of the system at 
time t and point x G M"^; by the analogy with the classical problem of the gas flow, 
it is natural to call these functions the density p{t,x) and the velocity v(t,x). A 
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nearly universal equation of continuity states that 

(1.4) npt + div(pi;) = a, 

where a = a{t, x, p) is the density of sources and sinks, and k is the fraction of the 
space available to the system. Moreover, the underlying physical principles dictate 
the equation of motion: 

(1.5) V = F(t, X, p, gradp) 

for the known vector- valued function F. An example is Darcy ' s Law for the ideal 
gas flow, 

(1.6) V = gradp, 

where p = pop"", a > 0, is the pressure, p is the permeability of the medium, and u is 
the viscosity of the gas. 

Assume that F{t, x, p, grad p) = — g(p)gradp — atplp) for some known functions q and 
ip and a constant vector a. Substituting (11.51) into (II. 4p . we then get 

(1.7) Kpt = A<l>(p) + a-grad^(p) + a(t,x,p), 

where $(x) = yq{y)dy and \['(x) = xip{x). Equation (11.71) appears in a variety 
of problems, including population dynamics and gas and water propagation. The 
underlying physical setting usually requires the solution p of (11. 7p to be non-negative 
for alH > and a; G M. 

The following two particular cases of (11.71) are worth mentioning: 

(1) Darcy's Law for the ideal gas flow (11.61) with constant p, z/, k and with a = 
results in equation (11.11) with m = 1 + a for the normalized (dimensionless) 
and suitably scaled density u of the gas. The main difference between (11.11) 
and the heat equation [m = 1) is that initial disturbance is propagated by 
(II. ip with flnite speed. 

(2) The same relation (II. 6p with constant p,^ and with k = const, a{t,x,p) = 
(To p, (To = const, leads to a basic model for the spread of a crowd-avoiding 
population: 

(1.8) ut = Au"" + aou, 

where m is a suitably normalized and scaled density of the population; see 
Gurtin and MacCamy |9,, Section 2]. The number ctq is the reproduction 
intensity of the population; the population is birth-dominant if do > 0. 

Comparing (11.20 and (II. 8p we conclude that (II. 2p can describe the dynamics of 
the crowd- avoiding population with random time- dependent reproduction intensity; 
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if X(t) = g{s)ds + f{s)dw{s), then 

(1.9) ao = g{t) + f{t)w{t), 

where w is Gaussian white noise. 

In general, allowing the function a in (11.71) to be random is a natural way to introduce 
randomness in the porous medium equation. If the positivity of the solution is not 
guaranteed with these values of a, the term Am™' in the equation is replaced by 
A(|-u|™~^-u). Similar to (11.21) . the randomness can be in time only, such as a{t, x, u) = 
+ uw in Mel'nik [131 [12] • The randomness can also be in time and space, such as 
a{t,x,u) = Ylik fk{'t)X)wk{t) in Kim [10], a(t,x,u) = F{u) + WQ(t,x) in Barbu et al. 
[2] and in Da Prato [5], or a(t,x,u) = F{u) + G{u)WQ{t,x) in Barbu et al. [3], with 
Wq representing Gaussian noise that is white in time but is sufficiently regular in 
space. Note that, unlike (II. 2p . none of the above models can benefit from Theorem 

o 

The derivation of (11.71) shows that another way to introduce randomness in the equa- 
tion is to allow velocity v to be random, for example, by considering a stochastic 
differential equation satisfied by v. Some results concerning the resulting porous 
medium equation were recently obtained by Sango fl6]. 

2. Nonlinear Equations with Multiplicative Noise 

If t> = v{t^x), t > X E satisfies the heat equation vt = Av, and c is a real 
number, then the function u{t,x) = v{t,x)e'^^ satisfies Ut = Am + cu. Similarly, if 
w = w{t) is a standard Brownian motion, then u{t,x) = v{t, x)e^^''^^~^^^'^^ satisfies the 
stochastic Ito equation du = Audt + udw{t). The objective of this section is to extend 
these results to some nonlinear equations. 

Consider the equation 

(2.1) Vt = F{v, Dv, D%, ...), t>0, X e 

with some initial condition. In (12.11) . v = v{t,x) is the unknown function, F is a give 
function, vt = dv/dt and D'^v denotes a generic k-th order partial derivative off with 
respect to x. 

We also consider the stochastic counterpart of (12.11) for the unknown random field 

u = u{t,x), t>0,xeR'^: 

(2.2) du = F{u, Du, D\ . . .)dt + u {f{t)dw{t) + g{t)dt), 

with the same initial condition as in (12. ip . In (12. 2p . / = /(/!:) is a locally square- 
integrable deterministic function, g is & locally integrable deterministic function, w 
is a standard Brownian motion on a stochastic basis F = (i7, JF, {JF^}j>q, P), and 
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the equation is in the Ito sense. We assume that F satisfies the usual conditions of 
completeness of JF and right continuity of Tt- 

Definition 2.1. Given a stopping time t, a classical solution of equation (12.21) 
on the set (0, r] = {{t,uj) : t < r} is a random field u = u{t,x) with the following 
properties: 

(1) u is continuous in {t,x) for all {t,uj) G (0,r] and x G M'^; 

(2) all the necessary partial derivatives of u with respect to x exist and are con- 
tinuous in {t, x) for all {t, uj) G (0, r] and x G W^; 

(3) The equality 



u{t, x) = u{0, x) + / F{u{s, x)^ Du{s, x), . . .)ds + / u{s, x){f{s)dw{s) + g{s)ds) 
Jo Jo 

for all {t, uj) G (0, r] and x G M'^. 

Taking in the above definition f = g = 0,we get the definition of the classical solution 
of equation (12. ip . It turns out that if the function F is homogeneous, then there is 
a one-to-one correspondence between the classical solutions of (12.11) and (12. 2p . The 
key component in this correspondence is a random time change. 

Definition 2.2. (a) We say that the function F is homogeneous of degree 7 > 1 

if, for every A > 0, 

(2.3) F{Xx,Xy,Xz,...) = X^F{x,y,z,...). 

(b) We say that equation (12. ip is homogeneous of degree 7 > 1 z/ the function F 
is homogeneous of degree 7. 

Proposition 2.3. Assume that the function F is homogeneous of degree 7. Define 

the functions 

(2.4) 

h{t) = exp (^j^ g{s)ds + f{s)dw{s) ~\ f\s)ds^ , H^{t) = h^-\s)ds. 

Then a function v = v{t, x) is a classical solution of (12. ip if and only if 

(2.5) u{t, x) = v{H^{t), x)h{t) 
is a classical solution of (12. 2p . 

Proof. Assume that f is a classical solution of (12. ip . Note that H'^if) = h''~^(t) and 

(2.6) dh{t) = h{t){f{t)dw{t) + g{t)dt). 
By the Ito formula, 

du{t,x) = Vt{H^{t), x)h'' {t)dt + u{t, x){f{t)dw{t) + g{t)dt). 
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Using (12.ip and homogeneity of F, 

Vt{H^{t),x)h^{t) = F(^h{t)v{H^{t),x),h{t)Dv{H^{t),x),h{t)D\{H^{t),x),...^ 
= F{u,Du,D\,. . .), 

and therefore m is a classical solution of fl2.2p . 

Conversely, assume that m is a classical solution of (12.21) . Since h{t) > for all 
t,LJ, the function is strictly increasing and has an inverse function R^. Define 
v{t,x) = z{R^(t))u{R^{t),x), where z{t) = l/h{t). Note that 

dz{t) = z{t) - g{t)dt - f{t)dw{t) + f{t)dt^ . 

Then, by the Ito formula, we conclude that f is a classical solution of (12. ip . □ 

As a simple illustration of Proposition 12.31 consider Burger's equation vt = vvx, t > 
with initial condition v{0,x) = —x. This equation is homogeneous of degree 2 and 
has a classical solution v{t,x) = —x/{l +t). Then the stochastic equation du = 
uUxdt + udw{t) with the same initial condition has a solution 

e"'(*)-(*/2) 
u(t,x) = —X -, . 

Proposition 12.31 can be generalized in several directions: 

(1) The functions f,g can be random and adapted. In fact, the process 
j^g{s)ds + f{s)dw{s) can be replaced with a semi-martingale X{t), possi- 
bly with jumps. Then h{t) becomes the stochastic (Dolean) exponential of X; 
see Liptser and Shiryaev [TTl Section 2.4]. 

(2) Other types of stochastic integral and other types of random perturbation 
can be considered, as long as the corresponding analog of equation ( 12. 6p can 
be solved and an analogue of the Ito formula applied. For example, consider 
the fractional Brownian motion and take f = 1, g = 0. With a suitable 
interpretation of the integral udB^ we get h{t) = (*)-(i/2)* ^ ggg^ for 
example, Nualart [TH Chapter 5]. 

(3) Since transformation (12.50 does not involve the space variable x, Proposition 
12.31 also works for initial-boundary value problems. 

(4) Transformation (12.50 can establish a one-to-one correspondence between gen- 
eralized (and even viscosity) solutions of (12. ip and (12. 2p . but the precise def- 
inition of the solution and the corresponding arguments in the proof require 
more information about the function F. 
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Let us emphasize that transformation (12. 5p does not lead to a closed-form solution of 
the stochastic equation fl2.2l) unless there is a closed-form solution of the deterministic 
equation (12. ip . Some methods of finding closed- form solutions of nonlinear equations 
of the type (12. ip are described in [6]. 

3. Stochastic Porous Medium Equation 

Recall that the classical porous medium equation is 

(3.1) = A(f/"^(t,x)), t>0, 

where Ut = dU/dt and A is the Laplace operator. This equation can model various 
physical phenomena for every m > 0; in what follows, we will consider only m > 1 
(m = 1 is the heat equation). We also assume that x G M'^ so that there are no 
boundary conditions. Note that, without special restrictions on m, the definition 
of the solution of (13. ip must include a certain non-negativity condition on U ; this 
condition is also consistent with the physical interpretation of the solution as a density 
of some matter. 

The scaled pressure V = V{t,x) corresponding to the porous medium equation 
(13.11) is defined by 

771 

(3.2) v(t,x) = ——U"'-\t,x) 

m — 1 

and satisfies 

(3.3) Vt = {m-l)VAV + \VV\\ 

where V is the gradient. The function V is extensively used in the study of the 
analytic properties of (13.11) . 

Let F = (fi, J^, {J-'t}t>o, P) be a stochastic basis with the usual assumptions, w = w{t), 
a standard Brownian motion on F, and r > 0, a stopping time. Let / = f{t) and 
g = g{t) be non-random functions such that / is locally square-integrable and g is 
locally integrable. 

Consider the following equation: 

(3.4) du{t, x) = A(M'^(t, x))dt + u{t, x){f{t)dw{t) + g{t)dt), t > 0, x G M"^. 

Definition 3.1. A non-negative, continuous random field u = u{t,x) is called a 
solution of equation (13.41) on the set (0,r] = {{tu}) : t < t} if, for every smooth 
compactly supported function ip = ip{x) the following equality holds for all (t,uj) G 
(0,rl.- 

(3.5) (M,<^)(t) = (M,vp)(0)+ / {u^,Aif){s)ds+ [ {u,ip){s){f{s)dw{s) + g{s)ds), 
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where 

{u,(p){t)= I u{t, x)(p{x)dx. 



Note that, with f{t) = g{t) = 0, this definition also apphes to the deterministic 
equation (13. ip . 

Define the functions 

h{t)=exp( [ g{s)ds+ [ f{s)dw{s)-l [ f{s)ds 
(3.6) ^ ^-^^ ^ 

H{t) = [ K^-\s)ds. 



Proposition 3.2. There is a one-to-one correspondence between the solutions of (13.11) 
and (13.41) given by 

(3.7) u{t,x) = U{H{t),x)h{t). 



Proof. Note that the porous medium equation (13. ip is homogeneous of degree m. 
Then (13. 7p is strongly suggested by Proposition [2l3l Since the solutions in question are 
not necessarily classical, the formal argument involves application of the Ito formula 
in the integral relation (13. 5p , but is completely analogous to the proof of Proposition 
[231 □ 



Two immediate consequences of (13. 7p are the comparison principle and maximum 
principle for equation (13. 4p : both follow from the corresponding results for the deter- 
ministic equation (13. ip . see the book by Vazquez [T7', Theorem 9.2]. 

Corollary 3.3. (a) Comparison principle: Ifu, u are two solutions of (13. 4p and 

ti(0, x) < n(0, x) for all x G W^, then u{t, x) < u{t, x) for all [t, u) G (0, r], x e R'^. 
(b) Maximum principle: If u is a solution of i\3A\) and < u{0,x) < M for all 
X G R"^, then < u{t,x) < Mh{t) for all {t,uj) G (0, r], x G M'^. 

Remark 3.4. 1. In the particular case f = 0, g = const, the relation (13. 7p was 
discovered by Gurtin and MacCamy [9J. 

2. While the current paper deals only with non-negative solutions of the porous 
medium equation, relation (13. 7p also holds for all the solutions of Ut = A.{\U\"^^^U) 
and du = A(\u\™'~^u)dt + udw{t). 

It is important to have the definition of the solution of (13. 4p on a random time 
interval, because even a classical solution of (13. ip can blow up in finite time, and 
then the corresponding solution of (13. 4p will blow up in random time. The quadratic 
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pressure solution provides an example. By direct computation, the function [/I^pI 
defined by 



J- I 12 \ 1 

\ m — I 



(3.8) f/['iri(t,x)= ,t„ = , g>0;ti = t 



is a classical solution of (13. ip for (t, x) e (0, tq) x M*^; [/^pI is known as the quadratic 
pressure solution because the corresponding pressure is 



(m- l)(t-t,)' 

for details, see Aronson [H pages 3-5] or Vazquez [TTl Section 4.5]. By Proposition 



«['iP](t,x) = f/[^Pl(i/(t),x)/i(t) 

is a classical solution of (13.41) on (0, r^J, where 

Tq = inf (t > : H{t) = tq) 

is the blow-up time; it is certainly possible to have < oo with positive probability. 
On the other hand, if / = 0, then, with a suitable choice of g, the life of every 
non-global solution of (13.11) can be extended indefinitely: note that g{t) = —a < 
corresponds to H(t) = (1 — e~^^)/a < l/a, and it remains to take a > sufficiently 
large. 

The following is the main result about global existence, uniqueness, and regularity of 
the solution of (13. 4p . 

Theorem 3.5. Assume that the initial condition m(0, x) is non-random and has the 
following properties: 

(1) non-negative and hounded: < u(0,x) <C,xE M'^; 

(2) continuous; 

(3) integrable: J u{0,x)dx = M, < M < oo; 

(4) square-integrable: J u'^{0,x)dx < oo. 

Then there exists a unique non-negative solution u = u{t,x) of (13. 4p . This solution 
is defined for allt>0,x& M.'^ and has the following properties: 

(1) the mean total mass satisfies E J u{t,x)dx = Me^^o aMd-^ ^ t > Q; 

(2) u{t, ■) is Holder continuous (as a function of x) for every t > 0; 
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(3) if {J^auP{0,x)dxy^^ = Mp<oo, then 
(3.9) (eJ uP{t,x)dx^ < Mpexp (^J^ g{s)ds + ^^^-^ f{s)ds^ . 

In addition, if the functions f, g are locally bounded, then the function u is Holder 
continuous on [T, +oo) x for every T > 0. 

Proof. For the deterministic equation, Sabinina fT5] proved existence, uniqueness, and 
conservation of mass (see also the book by Vazquez Chapter 9]; Caffarelh and 
Friedman [1] proved Holder continuity. It is also known [T71 Theorem 9.3] that 

UP(t, x)dx < / UP{0, x)dx, t > 0, p > 1, 

where U is the solution of (13.11) . It remains to use Proposition 13.21 and note that 

EhP{t) = exp (p [ g{s)ds + ~ [ f\s)ds^ , t > 0. 



Also, if the functions /, g are locally bounded, then the function h is Holder continuous 
of any order less than 1/2. 

□ 

Note that the initial condition of the quadratic pressure solution, 
g,i/(m-i) |^|2/(m-i)^ is not bouuded. Theorem 13.51 shows that a blow-up of the solution 
can be avoided with suitable growth restrictions on the initial condition. These re- 
strictions are sufficient but not necessary: consider, for example the linear pressure 
solution 

/m - 1 \ ^/^"""^^ 

u'^^\t,x) = max {H{t) + x, 0) h{t), t>0, xeR. 

By analogy with (13.21) . define the scaled pressure corresponding to equation (13.41) 
by 

771 

(3.10) v(t, x) = M"-i(t, x). 

m — 1 

An application of the Ito formula shows that v satisfies 

dv=( (m - l)v\v + \Vv\^ (m - l)(m - 2) ^ ,\ 

(3.11) V 2 'J 
+ (m — l)v {fdw + gdt). 
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On the other hand, equation fl3.3p is homogeneous of degree 2, and so Proposition 
suggests an alternative definition: 



(3.12) v{t,x) = V{H{t),x)h{t), 

where 



h{t) = exp (m - 1) / f{s)dw{s) - - — / f\s)ds 



2 jQ 



(m — l)(m — 2) I 2/ 



+ '"^ ^''"^ / f(s)ds + (m-l) / g{s)ds 
^ Jo Jo 

and H{t) = Jl^h(s)ds. An observation that ^(m — l)(m — 2) — (m — 1)^^/2 = 
— (m — l)/2 shows that, in fact, H{t) = H{t) and 

V(H(t),x)h(t) = U'^-\H(t),x)h"'-\t). 

m — 1 

In other words, f l3.10p and fl3.12p define the same function. Another way to phrase this 
conclusion is to say that, for the porous medium equation fl3.ip . the change of variables 
defined by Proposition [23] commutes with the transformation f/ i— >■ (m/ (m — !))[/''"" ^. 

4. BARENBLATT'S SOLUTION AND LONG-TIME BEHAVIOR 

Barenblatt's solution of Ut = A(f/™) is 

1 / / 1 I |2 \ \ 

(4.1) f/[BT](i,^.5) = _|'^ax(^6-^/?|L, ojj , t>0, XGM^ 
where b > and 

(m — l)c( + 2 

For the derivation of U^^^^ see Aronson PQ pages 3-4], Evans [3 Section 4.2.2], or 
Vazquez [T7[ Section 4.4.2]. The function has ^/P"""' the following properties: 

(1) the total mass of the solution, f U^^^^(t,x;b)dx, does not depend on t and 
is uniquely determined by b: by direct computation, 

(4.2) / [/PT](i^3;.5)rf^ = 5l/(2/3(n^-l)) U-l^ 



' Vm— 1 2/ 

where F is the Gamma function (see also Aronson [H Pages 3-4] and Vazquez 
[m Section 17.5]); 
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(2) t/P^l is a classical solution of Ut = A(f/"^) in the region 



(4.3) 

(3) For every p, q, > 0, Xq G M'^, the function 



(4.4) Uit, x) = i^^j f/PTl(pt + to, qx + xo; b) 

is also a solution of (13. ip . 
By Proposition 12.31 the function 

(4.5) u^''^\t,x;b) = U^''^\H{t),x;b)h{t), 

with H, h given by (13. 6p . is a solution of the stochastic porous medium equation (13. 4p . 
In particular, 

Barenblatt's solution U^^^^ determines the long-time behavior of every global solution 
of the deterministic equation (13.11) . Similarly, with obvious restrictions on H, uP"""! 
determines the long-time behavior of the solutions of the stochastic porous medium 
equation from Theorem 13.51 

Theorem 4.1. Assume that limt^infty H{t) = +oo with probability one. Let u = 
u{t,x) be the solution of the stochastic porous medium equation constructed under the 
assumptions of Theorem 1 5*. 51 Then, for every x G M*^, 

(4.6) lim {H{t)f''\u{t, x) - mP^] (t, x;b)\ = 



l/(m-l) 



t— >oo 



with probability one, where b is such that 

(4.7) / u{0, x)dx = 6V(2/3(-i)) ( !^ p\ ^ fe) 



2nm J r (^ + ^) ■ 

' \m—l 2/ 

Proof. This follows from Proposition 12.31 and the corresponding result for the deter- 
ministic equation, 

lim t^'^|f/(t,a;) -t/PT](t,x;6)| =0 

4— »oo 

(see Friedman and Kamin [S]). □ 
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Remark 4.2. Clearly, if lim^^oo H{t) is finite, then the long-time behavior of the 
solutions of fl3.4p will be quite different and, in particular, not as universal. By 
definition, the function H is non- decreasing and therefore has a limit (finite or 
infinite) with probability one. If f'^{s)ds < oo, then H(t) | +00 as long as 
limmft^oo Jq g{s)ds > —00. If f'^{t)dt = +00, then, by the law of iterated log- 
arithm, lim(_»oo H(t) can be finite with probability one, and more information about 
fit) and git) is necessary to proceed (see the example below). 

Further information about the asymptotic behavior of the solution can be obtained 
under additional assumptions about the functions /, g. 

Theorem 4.3. Assume that f'^{t)dt = 2a'^ and g{t)dt = fi for some a, /i G M. 
Then, for every x G M'^, 

(4.8) lim \u{t, x) - e« ^7^^! {e^"'-^^k, x;b)\=0 

>oo 

with probability one, where b satisfies (14.71) and is a Gaussian random variable with 
mean fi — and variance cr^ . 

Proof. Under the assumptions of the theorem, we have, with probability one, 

(4.9) \im (^J^ g{s)ds + f{s)dw{s) -^J^ f{s)ds^ = ^ 
Then lim^^oo h"^~^it) = e*^™~^^^ and therefore 

lim - / /i(™-^)(s)rfs = e('"-^)«. 

t^OO t Jq 

The result now follows by continuity of the function . Note that no rate can be 
specified in (14. 8 p without assumptions about the rate of convergence in (14.90 . 

□ 

Remark 4.4. It follows that equation (13. 4p DOES NOT have a nontrivial invariant 
measure in any traditional function space. Indeed, by (14. 8p . if limt^oo h{t) exists, 
then, for every x G M°', 

lim u{t, x) = lim f/^'^l (t, x; b) = 0, 

that is, either the solution decays uniformly or the mass spreads out to infinity. On the 
other hand, if limt^oo h{t) does not exist, then, by (13. 7p . no non-trivial limit ofu{t,x) 
can exist either. The same remark applies to the equation in a bounded domain. 

Other information about certain global solutions of the stochastic porous medium 
equation can be obtained by comparison with Barenblatt's solution. 
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Theorem 4.5. Assume that the initial condition u{0, x) of (13. 4p is continuous, non- 
negative, and compactly supported in M'^. Then, with probability one, 

(1) the solution u{t, x) is non-negative and has compact support in Mf^ for allt > 0; 

(2) the interface, that is, the boundary of the set {x e M*^ : u{t,x) > 0}, is 
moving with finite speed. 

(3) if limt^oo H{t) < oo, then the support of the solution remains bounded for all 
t > 0. 

Proof. By the maximum principle (Corollary 13.3( b) ) . if u{0, x) > 0, then so is u{t, x). 
Furthermore, if u{0,x) = for \x\ > R and u{0,x) < C, then, for sufficiently large 

Ci, 

u{0,x) < max(Ci — |x|^,0). 

By comparison principle (Corollary 13.3( a)) we then get u(t,x) < U{H{t),x)h{t), 
where U{t,x) is a function of the type (14. 4p with p = q = to = 1, xq = 0, and 
b sufficiently large. Therefore, u{t,x) = for |x| > C2H^{t) for a suitable (non- 
random) number C2. □ 

Example. Consider the equation 

(4.10) du = Au'^dt-^udw{t) 

and assume that m(0, x) is continuous, non-negative, compactly supported, and 
J^au{0,x)dx > 0. Then there exists a random variable rj such that < rj < 00 
with probability one and u{t, x) = for all |x| > rj and all t > 0. Indeed, in this case 

/i(t) =e'"W-(*/2), H{t)= f e^^'^-^'/^^ds, 

Jo 

and, by the previous theorem, it is enough to show that H{t) is bounded with prob- 
ability one. Let Tq be the last time w(t) exceeds t/4: Tq = sup{t > : w{t) > it/4}. 
Then 

Hit) < r e"'«-(*/2)rft + 4e'ro/4^ ^ > g. 
Jo 

By the law of iterated logarithm, Tq < 00 with probability one, and therefore 
limt^oo H{t) < 00 with probability one. 

Notice also that lim^^oo h{t) = with probability one, and consequently, for every 

X e M"', 

lim u{t,x) = lim h{t)U (H (t) , x) = 0, 

t—*oo t— >oo 

because the solution U = U{t, x) of the deterministic equation (13. ip with initial 
condition U{0, x) = u{0, x) is a uniformly bounded function. On the other hand, 
we know that J^j^U(t,x)dx = J^au{0,x)dx, and, since ¥,h(t) = 1 for all t > 0, 
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we conclude that E J^j u(t,x)dx = f^au{0,x)dx. In other words, the solution of the 
stochastic porous medium equation fl4.10p is supported in the same (random) compact 
set for alH > and decays to zero as t — > oo, while preserving the mean total mass. 
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